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Abstract 

We discuss the relation between singularities of correlation functions and causal 
properties of the bulk spacetime in the context of the AdS/CFT correspondence. In 
particular, we argue that the boundary field theory correlation functions are singu- 
lar when the insertion points are connected causally by a bulk null geodesic. This 
implies the existence of "bulk-cone singularities" in boundary theory correlation func- 
tions which lie inside the boundary light-cone. We exhibit the pattern of singularities 
in various asymptotically AdS spacetimes and argue that this pattern can be used 
to probe the bulk geometry. We apply this correspondence to the specific case of 
shell collapse in AdS/CFT and indicate a sharp feature in the boundary observables 
corresponding to black hole event horizon formation. 
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1 Introduction 



The AdS/CFT correspondence, which relates quantum gravity in asymptotically Anti de Sit- 
ter spacetimes to a non-gravitational gauge theory, is an invaluable window to understanding 
effects of gravity beyond the semi-classical approximation. In principle the correspondence 
provides a framework to address many long standing questions in quantum gravity, such 
as the nature of singularity resolution or issues relating to the information paradox. While 
much has been learnt in the past few years, the status of the dictionary between gravita- 
tional quantities and field theory observables is still at a somewhat rudimentary stage. This 
has been a stumbling block for exploiting the correspondence to its full power, especially in 
eliciting answers to oft posed quantum gravitational questions. 

In classical general relativity, concepts such as causal structure, event horizons, singular- 
ities, etc., play an important role in understanding the geometry of the spacetime manifold. 
In the semi-classical approximation these concepts are useful in understanding the dynamics 
of quantum fields in curved backgrounds. Given the AdS/CFT correspondence it is inter- 
esting to understand field theoretic encoding of these geometric concepts. Considering the 
central role played by geometry in classical general relativity, one naively expects them to 
have a well defined representation in the field theory. 

For instance, consider the field theoretic representation of bulk causal structure. In the 
semi-classical limit, the bulk causal structure can be read off from the properties of correlation 
functions of quantum fields propagating on the spacetime manifold. For asymptotically AdS 
spacetimes the bulk correlation functions in turn determine (in some suitable scaling limit) 
the boundary correlation functions. Consistency within the correspondence requires therefore 
that the boundary correlation functions are compatible with the bulk causality constraints. 
This has been demonstrated to be true in many examples [1-4]. A crucial ingredient in 
most of these analyses has been to exploit the asymptotic AdS geometry where consistency 
is guaranteed by the symmetries in question. 

The non-trivial aspect of the geometry that one would like to understand is the behaviour 
of the bulk causal structure deep in the interior of an asymptotically AdS spacetime. For 
example one can ask whether the presence of an event horizon in the spacetime has a non- 
trivial signature in the boundary. The zeroth order answer to this question, at least for 
eternal black holes in AdS spacetime, is that the dual field theory is in a thermal state. 1 So 
the presence of the horizon is encoded in a new scale for physics in the boundary theory, the 
thermal scale. However, this information hardly probes the causal structure in the vicinity 
of the horizon. 

Furthermore, since the AdS/CFT correspondence is most firmly understood in the Eu- 
clidean framework, questions pertaining to genuinely time-dependent processes in the bulk 

1 This is true only for black holes whose horizon size is larger than the AdS scale. 
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and holographic duals thereof become much more subtle. For example, consider modelling a 
time dependent process, such as of black hole formation in the bulk, within the field theory. 
A more nuanced question, then, would to ask whether there is any field theory information 
to be gained about the event horizon formation process. Phrased differently: can we see 
horizon formation directly in the gauge theory? Since in the bulk the event horizon forms at 
a sharply-localized event 2 (despite the fact that we need to know the full future evolution of 
the spacetime to actually find the horizon), we would expect that this will manifest itself in 
some correspondingly sharp feature in the gauge theory. 

The conventional lore is that the UV/IR relation, which maps local regions in the interior 
of the bulk to non-local objects in the CFT, makes it hard to extract useful information about 
bulk geometry (c/., [5] for attempts to use UV/IR relation to extract information about the 
horizon). Given that the geometric scale associated with the horizon corresponds to the 
thermal scale for neutral black holes, it is hard to see how to extract a precise signal about 
the causal structure in the neighbourhood of the event horizon. 

The clue comes from the progress made in the recent past in identifying the CFT signature 
of the black hole singularity [6,7]. In this case the field theoretic observables transcend the 
classical barrier of the event horizon to encode information of behind the horizon physics. 
The main observation of [6], which built on the original work of [8,9] was to use the intrinsic 
non-locality of the boundary correlation functions to identify signals of the singularity. The 
black hole singularity was identified to correspond to a particular 'light-cone' like singularity 
in the field theory correlation function. 3 It has further been proposed that this technique 
can be exploited to understand physics behind the event horizon and in particular used to 
investigate aspects of inflationary geometry within AdS/CFT [10]. The basic idea was to 
look at the bulk Green's functions in a saddle approximation where they are dominated 
by geodesies and use the bulk computations to derive predictions for boundary correlation 
functions. 

In the present work we argue that the horizon formation in gravitational collapse can 
be detected in the boundary theory by examining the structure of singularities of generic 
Lorentzian correlators. In particular, we show that a sharp horizon-formation time can be 
extracted from the pattern of singularities. 

Our argument is based on a connection between null geodesies in the bulk spacetime and 
singularities in correlation functions of local operators evaluated in the state corresponding to 

2 In ths work we consider spherically symmetric spacetimes. More generally, in non-spherical spacetimes 
the event horizon could form at a locus of events; our methods should generalise to these cases as well. 

3 The identification of the singularity in the strict large N limit is easiest when formulated in terms of 
momentum space correlators [7] as opposed to direct computation in position space where the singularity is 
not visible in the primary sheet of the correlator [6]. In particular, it was confirmed in [7] that the signatures 
of the singularities go away when the rank of the boundary theory gauge group is finite, implying that the 
singularities are resolved in quantum gravity. 
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the bulk geometry. 4 The basic idea is that CFT correlators will exhibit light-cone singularities 
when the points of operator insertion are connected by a null geodesic. The connection 
implies that CFT correlators in excited states have additional Lorentzian singularities inside 
the light cone, which we will call bulk-cone singularities. Usually CFT correlators exhibit 
light-cone singularities and the location of these is determined by the causal structure of 
the background spacetime on which the field theory lives. For a CFT living on R x S d_1 
these would be given by the conventional light cone of the Einstein static universe. Our 
main observation is that bulk causality, together with the dictionary between bulk Green's 
function and boundary correlation functions in the AdS/CFT context, necessitate additional 
bulk-cone singularities in the CFT correlators. 

Armed with this relation we can ask how certain geometric structures are encoded in 
the boundary data. In addition to a black hole collapse geometry which motivated this 
study, we also look at the geometry of a star in AdS spacetime and the eternal black hole 
geometry. In all these cases, interesting details about the geometry and its causal structure 
that can be read off from the properties of the singularities in the CFT correlators. In fact, 
as has recently been shown in [11], the location of the boundary singularities can be used to 
reconstruct the complete bulk metric for a class of static, spherically symmetric spacetimes. 
Our analysis is carried out by explicitly studying geodesies in the bulk spacetime and using 
this to infer properties about correlation functions in strongly coupled gauge theories. While 
there are a few consistency checks we can establish in some simple cases, many of our results 
are predictions for strong coupling CFT correlation functions in certain excited states. 

The plan of the paper is as follows: We begin in Section 2 by motivating our claim that 
Lorentzian AdS/CFT correspondence implies that correlation functions of local operators 
evaluated in excited states of the field theory exhibit light-cone singularities whenever the 
operator insertion points are connected by null geodesies through the bulk spacetime. We 
then confirm the relation between light-cone singularities and bulk null geodesies for vacuum 
correlators. Further, we generalize to excitations about vacuum state and obtain predictions 
for the singularities in the correlation functions from the behaviour of null geodesies. To 
demonstrate how sensitively the singularity pattern depends on the state, in Sections 3 and 

4 we explicitly study the nature of the singularities and their implications for correlation 
functions in two scenarios: a generic excited state and a thermal state, which we model in 
the bulk by a star in AdS and an eternal black hole geometry, respectively. We devote Section 

5 to our main example, the collapse of a black hole in AdS spacetime, and elicit from the 
geometry signatures of event horizon formation in the field theory. We conclude in Section 6 
with a discussion. Finally, in the Appendices we present more detailed calculations pertaining 
to the various geometries considered (pure AdS, star in AdS, eternal Schwarzschild-AdS, 
Poincare patch, null shell collapse, and Vaidya-AdS). 

4 Spacetimes which are asymptotically AdS can be thought of as deformations of pure AdS by normalizablc 
modes in the supergravity description. In field theory terms this corresponds to a generic excited state 
obtained by acting on the CFT vacuum with the appropriate operator. 
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2 Null geodesies and boundary singularities 

In this section we motivate a simple relation between singularities of the boundary correlation 
functions and null geodesies in the bulk spacetime. The main idea is to exploit the AdS / CFT 
dictionary for evaluating strong coupling correlation functions in the CFT using bulk Green's 
functions. Since the bulk Green's functions are sensitive to the bulk causal structure we 
expect that the properties of the bulk light cone to be directly visible in the boundary 
correlation functions. 

A short word on notation: for convenience, we will denote the d + 1 dimensional bulk 
spacetime by M and its d dimensional boundary by DM. Further, we use labels x,y etc., 
for points on the boundary dM and r to denote the radial coordinate normal to dM. 



2.1 General argument 

The AdS / CFT correpsondence gives us a precise map relating string theory (or its low energy 
supergravity limit) on M with a field theory living on the boundary dM of M. To be precise, 
in the low energy supergravity approximation of the bulk theory (which is good as long as 
we choose the 't Hooft coupling A and the rank of the gauge group N on the boundary to 
be large, i.e., A ^> 1 and iV ^> 1) we can obtain the boundary correlation function G(x,x') 
from the bulk propagator Q(x, r; x'; r') of the corresponding field by taking both end points 
to the boundary. More precisely, 

G(x, x') — 2 v lim (rr') A Q(x, r; x'\ r') , (2-1) 

r— >oo,r'— *oo 

where Q is the free field theory propagator for the corresponding bulk field. Focussing on 
the simple case of local gauge invariant operators 0(x) of conformal dimension A in dM 
which are dual to free scalar fields parameterized by mass m in the bulk, in (2.1) we have 

A = ^ + z/, u= ^m2 + ? , (2.2) 

where we have set the curvature radius of AdS to be 1. Note that by choosing the asymptotic 
AdS time to coincide with the boundary time, one ensures that time ordering chosen for the 
bulk propagator Q carries over to the boundary correlator in the limiting procedure. In what 
follows we will consider Feynmann propagators. 

That the limit in (2.1) is well-defined for an asymptotically AdS spacetime can be seen 
as follows. Upon fixing (x', r'), Q(x, r; x', r') is a normalizable solution of the free bulk wave 
equation on M. This implies for an asymptotically AdS spacetime 

Q{x, r; x', r') ~ r~ A , r — * oo . (2.3) 

Similarly, 

G(x, r; x', r') ~ r'~ A , r' — > oo . (2.4) 
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Thus for generic boundary points x and x', the limit (2.1) is well defined and yields a regular 
G(x,x'). However, if x and x' are connected by a null geodesic, (2.3) and (2.4) could break 
down, in which case one expects 

lim Q(x, r; x'\ r') -U . (2-5) 

r,r'—* oo 

A naive application of (2.1) will yield a divergent answer as the limit is not well-defined. 

In order to obtain explicitly the singular behavior of G(x, x') has to work with a suit- 
ably regularized expression. To this end one can first consider some x" lying in a small 
neighborhood of x 1 in the boundary. If x" is not connected to x by null geodesies, (2.1) can 
be straightforwardly applied to obtain G(x,x"). One can then take x" — > x' to obtain the 
desired correlator G(x, x'). 

We will now motivate the above discussion in a slightly different perspective using the 
geodesic approximation to Q(x,r; x' ,r'). For this purpose, we will take m large 5 so that 
A ~ m. To consider the limit (2.1), we put the boundary dAi at a cut-off surface r = 
r' = A — * oo, on which (for an asymptotically AdS geometry) the induced metric can be 
approximated as 

ds 2 = A 2 ds 2 bd (2.6) 
where ds 2 d denotes the metric on the boundary dAi. 

Consider points A = (x, A) and B = (x', A) on the cut-off surface, which are spacelike 
separated. In the large m limit, one expects that the saddle point approximation is valid. 
Denoting by d(A, B) the proper distance between two (spacelike separated) points in Ai we 
have 

Q(A, B) oc e~ mdiAB) ~ e~ 2m logA , (2.7) 

where we used the fact that as A — > oo, the proper distance d(A,B) ~ 2 log A + • ■ ■ in 
asymptotically AdS spacetimes. Note that (2.7) is consistent with (2.3) and (2.4) and hence 
makes sure that the limit (2.1) is well defined. 6 

Now suppose that there is a null geodesic connecting A and B - in this case d(A, B) = 
and Q(A, B) ~ 0(1). As suggested above we consider points C = (x",A) in a small 
neighborhood of B, for which 

d{A,C) ~log[A 2 5x 2 ] (2.8) 

where 5x is the proper distance between x and x" on the boundary (i.e., in terms of metric 
ds 2 d in (2.6) on dAi). Plugging (2.8) into (2.7) and then taking the limit as suggested in 
(2.1), we have 

G(x,x")~-^- 7 z"->z\ (2.9) 

5 The relation between null geodesies and singularities of G(x,x') does not require taking m large. 

6 For generic timelike separated points A,B on dA4, we also expect (2.7) to be true based on analytic 
continuation from spacelike separated points. Note that in the pure AdS, Q is a function of the proper 
distance d(A, B) which can be analytically continued in the complex plane to timelike separation. 
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To summarize, we have argued that a boundary correlation function G(x, x') = (O(x)O(x')) 
is singular if and only if 7 there exist null geodesies connecting x and x' . Our argument is not 
completely rigorous and we discuss potential loopholes below. This proposal can be checked 
in simple examples in which one can work out the boundary G(x, x') explicitly, like global 
pure AdS, Poincare patch and the BTZ black holes (see appendices). 



2.2 Types of null geodesies and causal structure 

In this subsection we discuss the location of the singularities in G(x, x') and structure of null 
geodesies. For this purpose, let us fix x and consider those x' which are connected to x by 
null geodesies. The null geodesies can be divided into two types: those lying entirely in the 
boundary (type A) and those lying in the bulk except for their end points (type B). Type 
A geodesies are simply null geodesies of the boundary theory and define the light cone of 
the CFT. For points x and x' connected by Type A null geodesies, G(x, x') has the standard 
light-cone singularities. Type B null geodesies are more interesting and we now discuss them 
in some detail. 

In pure global AdS (where the boundary is R x S d_1 ), all Type B null geodesies connect 
the anti-podal points on S d ~ l . The points in question are given by (ti, Q) and (t Q , — Q), with 
t — ti = 7r R-Ads — 71 an d ±^ denote points which are anti-podal on the S d ~ l . Note that 
the (boundary) time separation t a — U = ir needed for the geodesic to traverse the bulk is 
exactly the same as that of a Type A null geodesic connecting the anti-podal points only in 
pure AdS spacetime; the former becomes longer for general asymptotically AdS spacetimes. 
Further details on geodesies in AdS, both analytic and pictorial, appear in Appendix A. 

Since in the case of pure AdS the end points of the Type B null geodesies already coincide 
with the boundary light cone, they do not give rise to additional singularities in the boundary 
correlation functions. It is easy to verify that the behaviour of null geodesies in AdS is 
consistent with the singularity properties of boundary correlation functions. The boundary 
correlators in questions are two-point functions in the CFT vacuum of primary operators of 
dimension A which for a d-dimensional CFT on R x S d ~ l are given by 

(O (t, fi) O (t', fi')> oc - 1 , (2.10) 

(cos(t — r) — cos(iz — U')) 

The singularities of the correlation function (2.10) are given by 5t = 5Q with 5t = t — t' , 
etc., which coincide with the light cone of the boundary manifold. 

In a perturbed AdS spacetime, the end points of Type B geodesies in general do not 
coincide with the light cone of the boundary manifold. They will therefore give rise to 



7 While we have explicitly shown that bulk null geodesies lead to new singularities in boundary correlation 
functions, our arguments indicate that the converse should also be true. This is equivalent to the statement 
that (2.5) is true only for points connected by null geodesies. 
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additional singularities in the boundary correlation functions. We will call such singularities 
bulk-cone singularities. An interesting feature of asymptotically AdS spacetimes is that 
there exists a lower bound on the 'time delay' experienced by geodesies exploring the bulk 
spacetime. In a situation where the bulk manifold A4 with a timelike conformal boundary 
dM. satisfies the null energy condition, the null generic condition, and strong causality 
and compactness, it has been proven by Gao and Wald [3] that the end points x' of Type B 
geodesies always lie inside the light cone of the boundary manifold, i.e., x' and x are timelike 
separated. This implies that any generic perturbation of AdS (satisfying the conditions 
mentioned) will produce a time delay of null geodesies relative to pure AdS. For example, 
imagine a state corresponding to a thermal gas in AdS at some internal mass density p Q . 
Then a radial geodesic will arrive to the anti-podal point of the boundary after a time At 
which increases monotonically with increasing p a . 

This implies that correlation functions in some excited states 8 of the theory will have 
additional Lorentzian singularities inside the light cone. We also note that it can be checked 
explicitly that for backgrounds with non-compact boundary, Type B null geodesies do not 
exist for states preserving the boundary symmetries (for a proof of the statement see Ap- 
pendix D). Thus additional singularities appear to arise only for CFTs on compact spaces 
or those that violate the Poincare symmetery for CFT on R 3,1 . 

While the theorem of Gao and Wald only applies to asymptotically AdS spacetime sat- 
isfying the conditions stated above, from the causality of the boundary theory one would 
like to conclude that the end points of Type B null geodesies should always lie inside the 
light cone of the boundary for any asymptotically AdS background with a field theory dual. 
Otherwise it would imply that one can set up a special state in the field theory to transmit 
signals faster than the speed of light (see also [12]). 

To conclude this section let us comment on a potential loophole in our argument. Our 
connection between null geodesies and singularities of boundary correlation functions is based 
on the validity of (2.5) for any two points connected by a null geodesic. While this assumption 
is consistent with our general understanding of quantum field theory in a curved spacetime, 
we do not have a rigorous proof of the statement. 9 Also note that one should distinguish be- 

8 In order to have non-negligible back reaction on the bulk metric, the ADM mass M of the perturbation 
should satisfies GnM ~ 0(1). This implies that the corresponding excited state in boundary theory should 
have energy of order 0(N 2 ). In the following discussion, by generic excited states, we mean generic states 
of such energies. 

9 For example, in a semi-classical approximation to the propagator, one may consider the following sce- 
nario: In general, the end points of a null geodesies can be shared by other spacclikc and/or timclike 
geodesies. If the propagator is given by summing over all geodesies, then the null geodesic clearly dominates. 
However, there might be a situation that path integration contours for the propagator cannot be deformed 
to the saddle point corresponding to the null geodesic. Then the null geodesic will not contribute to the 
propagator, even though it naively dominates. Analogous phenomenon was encountered in [6] for a spacelike 
(almost null) geodesic 'bouncing off' a black hole singularity; though it is not clear whether this contingency 
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tween points which can be connected by a null geodesic and those points which the geodesies 
connecting them only approach null. In the latter case, when there does not exist a null 
geodesic which connects two points, our argument does not directly apply. We will elaborate 
more on this when discussing the gravitational collapse background in Section 5. 



3 Static asymptotically AdS geometries 



We now turn to applying the results of Section 2 to static spacetimes. A trivial example 
where our relations can be verified is of course the empty AdS spacetime which we have 
already commented on. We now turn to consider geometries which are deformations of AdS, 
with a free parameter such that pure AdS is retrieved in a particular limit of this parameter. 
We will first consider situations with sufficient symmetries to simplify the analysis before 
turning to more complicated examples. For concreteness, we will work in 5 dimensions, i.e., 
take the spacetime to be asymptotically AdS$ (times a S 5 which will not play any role). 



3.1 Geometries of stars in AdS 

As our first nontrivial example, we will take our spacetime to be static and spherically 
symmetric, and to make things even simpler, we will take the matter content to be that 
of a self-gravitating gas of radiation. We note that self-gravitating radiation in AdS has 
been explored previously in [13] and related geometries of boson stars in asymptotically AdS 
spacetimes have been considered in e.g. [14,15]. The symmetries imply the corresponding 
stress tensor is that of a perfect fluid; radiation equation of state which makes the stress 
tensor traceless, i.e., 

T ab = p(r)u a u b + P(r) (g ab + u a u b ) , P = - p (3.1) 

where u a is a co-moving gas 4- velocity, and p(r), P(r) are density and pressure. Since the 
radiation will be confined by the AdS potential, we will refer to this configuration as a 
"star" in AdS. The metric can be obtained by solving the Einstein's equations with negative 
comological constant and the requisite stress tensor, as discussed in detail in Appendix B. 
Here we summarize the main results. 

The metric can be written in the form: 

ds 2 = -f{r)dt 2 + h{r)dr 2 + r 2 dfl 2 3 (3.2) 

with 



h(r) 



9 m\r) 
r 2 + 1 - 



(3.3) 



will ever arise for strictly null geodesies. Given that one can send light signals following null geodesies, one 
might be able to argue a null geodesic should always contribute to a propagator. 
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and 



f(r) 



Pc 



p(r] 



2/5 



where the mass function m(r) is defined in terms of an integral of the density p(r), 

, , 2 rr 



m r 



p(r) r 



(3.4) 



(3.5) 



and poo is the coefficient of the leading fall-off of p(r), determined by p(r) ~ Poo/r 5 as r — > oo. 
The field equations specify the system of coupled first order ODEs which m(r) and p(r) must 
satisfy: 

m'(r) = | p(r) r 3 (3.6) 

(3.7) 

with boundary conditions m(0) = and p(0) = po, where po is a free parameter of the con- 
figuration, specifying the internal density of the gas. Thus the geometries are parameterized 
by a single parameter p , with pure AdS retrieved in the po = case. 



P \ r ) ~ — — — m(r) — 




(a) 



m(r) 




2 3 



2 3 



Fig. 1: Density and mass functions for the "star" geometry, where the central density is set to po = 10. 



From (3.6) and (3.7) we find the following asymptotic behavior for m(r) and p(r) 



r — > oo : 



: 



( \ P™(\ 5 | 5 ( 4M + 7 ) 



m(r) — > M 



p(r) po + 0(r 2 ), m(r) -> ^r 4 + 0(r e ) 

o 



(3.8) 

(3.9) 
(3.10) 



where p^ and M (M is proportional to ADM mass) are constants which are determined by 
Po. Using (3.3) and (3.4) we also find that 

r — > oo : 



/(r)^ + l-- 



r 2 + 1 



r -> : 



/(0 



Poo 

Po 



+ ■ • ■ h(r) - 
+ 0(r 2 ), /i(r) -> 1 + 0(r 2 ) 
9 



(3.11) 
(3.12) 



One can readily solve the full equations (3.6) and (3.7) numerically. Fig. 1 shows the 
density and mass functions for the internal density po = 10. 



f(r) (a) h(r) (b) 




12345 12345 



Fig. 2: Metric functions for the "star" geometry with central density po = 10 (red curves), and for 
comparison corresponding metric functions in the pure AdS geometry (blue curves; f(r) is the 
higher and h(r) is the lower of the two curves). 

The metric functions are plotted in Fig. 2, again for the same set-up as in Fig. 1. For 
comparison we also plot the corresponding metric coefficients for pure AdS; as we would 
expect, the metric does approach that of AdS as the density becomes small. 

3.2 Geodesies in the star AdS geometry 

We will now turn to the properties of bulk null geodesies (Type B) in the star geometry. 
Staticity and spherical symmetry implies that E (energy) and J (angular momentum) are 
conserved along any geodesic and hence we will parameterize the geodesies by these quanti- 
ties. The equations for null geodesies are: 

a 1 

i= — , and f(r)h(r)r 2 = a 2 -V(r) , (3.13) 

j(r) r z 

where 

V(r) = ^, a = | (3.14) 

and the dot denotes differentiation with respect to the affine parameter along the geodesic. 
For null geodesies the absence of scale implies that the relevant parameter is the ratio a = 
E/J. Useful quantifiers for figuring out which points on the boundary are connected by 
null geodesies are the temporal and angular separation of the geodesic endpoints, which we 
denote by At and Atp, respectively. Since we have a numerical solution for the metric, we 
have to integrate the geodesic equations numerically to find the desired properties. 

Consider first radial null geodesies, which we expect will emerge at a later time then the 
corresponding geodesic in pure AdS, consistent with the theorem of [3]. Denoting this time 
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by Ato, we have 



At 



o 




h(r) 



dr 



<o i, f(r) 

For a star geometry with central density po = 10, we find numerically that At 

At 



(3.15) 
3.696, 



At 

1 . 3 
1 .25 

1.2 
1 .15 



AdS 



1 .05 



20 40 60 



:0 100 



pO 



Fig. 3: Time delay for radial null geodesies through star in AdS, as a function of the star's internal 
density p$. 



which is 1.177 times longer than in pure AdS. Fig. 3 exhibits the dependence of this time 
delay on the value of po- We see that while the time delay increases monotonically with 
increasing density, it does so more and more slowly. 10 

Now let us consider null geodesies with nonzero angular momentum. (Appendix B further 
explores non-radial spacelike geodesies.) In Fig. 4, we show null geodesies projected onto a 
spatial (constant t) slice (left) and the t — r plane (right). Not surprisingly, the star has a 
focusing effect on the geodesies. The important point to note is that now, unlike the pure AdS 
case, there is a finite spread of the At and A<p endpoint values the null geodesies can take. 
More specifically, because of the attractive nature of the star's gravitational potential, the 
small angular momentum (J/E <C 1) geodesies 'overshoot' in terms of angular separation, 
and at the same time they experience some time delay, which means that both At and A<p 
increase. However, this effect does not keep increasing monotonically with increasing angular 
momentum, because the further from the origin the geodesies penetrate, the less effect of 
the star they feel. In the high angular momentum limit (J/E ~ 1), the geodesies hug the 
boundary and consequently behave just as the corresponding geodesies in pure AdS. 



10 The intuitive reason for Ato not increasing more rapidly with po is that the larger the internal density 
po, the faster the density function p(r) falls off, in such a way as to keep the total mass M bounded, as 
discussed in Appendix B. Numerically we find that Ato appears to be bounded as po — > oo. 
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Fig. 4: Null geodesies in star with po = 10 in AdS, projected onto a constant t slice and the t — r 
plane, for varying angular momentum to energy ratio (E = 10 and J = 0, 1, ... , 10). On the 
left, the bold circle corresponds to the AdS boundary, whereas on the right, the dashed vertical 
line corresponds to the origin r = and the bold vertical line to the boundary at tanr = ^. 
The range of t plotted is (0, l.lAirj)- (Analogous plots with larger po appear in Fig. 19 of 
Appendix B.) 

To illustrate explicitly how the endpoints of these various geodesies compare, Fig. 5 sum- 
marizes the behaviour of At and Atp for the various geodesies. In (a), where the endpoints 
are plotted for small internal densities, po = 1, . . . , 10, we see that the effect of overshooting 
becomes more pronounced the higher p is. In (b), where po = 10, ... , 100, we illustrate that 
the effect however remains bounded even for high internal densities. (We have verified this 
numerically for po ^ 10 18 .) Note that the At intercept for each curve can be read off already 
from Fig. 3. 

The boundedness of At implies one of the interesting features for the star geometries, 
namely the absence of null geodesies orbiting the star. In fact, one can obtain this result 
more directly, as follows: From (3.13) we see that in order to have null geodesies at fixed 
radial distance from the star, i.e., r(A) = constant, we require that a 2 — V(r) = and 
4-V(r) = 0. It is easy to show numerically that these conditions are never satisfied for the 
star geometries. 11 As we will argue below, this absence of circular orbits around the star 
leads to marked differences in the behaviour of singularities for correlation functions in the 
pure state corresponding to the star and in themal density matrix represented by an eternal 
black hole geometry. In particular, in the star case, the spread of the geodesic endpoints 
remains bounded, unlike for the black hole case which has a circular orbit. 

11 The two conditions can be distilled into a necessary condition for circular photon orbits; defining Q(r) = 
j2 p{t) r 2 + 2 ^jP- we need Q(tq) = 1 at some r = r^. We find that Q(r) < 0.6 for all r irrespective of the 
internal density pg. 
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Fig. 5: Endpoints of null geodesies in the AdS geometry of star with (a) po = 1,2,..., 10 and 
(b) po = 10,20, . . . , 100. Each curve is plotted by varying a = -j. The top of the curves 
corresponds to big values of a; in fact, one can easily show that the slope of each curve (as a 
function of a) is given by I /a. 

The essential points relating to null geodesies in the star background can be summarized 
as follows: 

• The radial null geodesies which go through the bulk geometry emerge at the anti-podal 
point on the sphere (A<p = ir) at a later time than in pure AdS (At > 7r). The time 
delay is a monotonic function of the internal density po- 

• Null geodesies with angular momentum exhibit both a time delay as well as a shift in 
Aip. For small angular momentum, both At and Aip increase with increasing angular 
momentum, whereas for large angular momentum they decrease to match with the 
AdS value At = A<p = n when J = E. While this effect increases with increasing p , 
At and Aip remain bounded for arbitrary po- 

• There are no null geodesies at fixed value of the radial coordinate i. e., the star spacetime 
does not admit circular photon orbits. 

We stress that since the endpoints of the null geodesies indicate the location of a singu- 
larity in the corresponding correlation function as discussed above, the data of Fig. 5 is in 
principle easily extractible from the gauge theory. This demonstrates how from the gauge 
theory we can read off the details of the geometry deep inside the bulk. In other words, just 
by looking at the spread of the endpoints of the null geodesies, we can distinguish the star 
from a small black hole (as demonstrated below) or other configurations of the same mass 
in AdS. Note that this is rather different method of extraction of details of deep IR in the 
bulk from that employed in [16]: the latter considered one-point functions, whereas here we 
make explicit use of the nonlocal nature of the two-point functions. 
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3.3 Field theory signature of excited AdS geometries 

We now turn to a brief analysis of the star geometry from a field theory perspective. By 
construction we have an asymptotically AdS geometry which is static and spherically sym- 
metric. While we have modeled the star by an effective equation of state, it is easy to see 
that similar solutions can be obtained in gauged supergravity theories. In any event the 
geometries in question correspond to normalizable deformations of the AdS vacuum, im- 
plying that they are dual to states in the field theory. In particular, these states should 
be eigenstates of the CFT Hamiltonian (since they are static). One can therefore use the 
state operator correspondence in the CFT to map the state into an operator insertion; for 
convenience we will call the state of the star \star) and the corresponding operator O* i.e., 
| star) = |0). Note that to obtain a geometry with spherical symmetry we would need to 
smear the operator insertion 0* on the boundary S 3 making it effectively non-local. 

One of the interesting questions that we will be unable to provide a concrete answer for 
is the precise description of the state \star) in the field theory. From general arguments it is 
clear that the operator has dimension of 0(N 2 ) in the field theory. 12 In order to ascertain 
the precise state we need more information than the one point function of the stress tensor 
(which is of course given by the leading fall-off in the metric). However, it is clear that on 
physical grounds such states ought to be constructible in the field theory and we will proceed 
with the assumption that is a bona fide operator in the field theory. 

Our analysis of geodesies in the star spacetime can be rephrased in light of the discussion 
in Section 2 as predictions for new 'light-cone' singularities in the boundary correlators in the 
state | star). The later arrival at the boundary of radial null geodesies consistent indicates 
that correlation functions of generic operators O in the field theory, {star\ 0(x) 0(y) \star) 
will have singularities at Aip = it and At given by the appropriate time delay. Similarly, the 
implication of bulk Type B geodesies with angular momentum is that we will have light-cone 
singularities in field theory correlators at operator insertion points x, y such that t x —t y = At 
and Q x — Q y = Aip with At and Aip given by the geodesic analysis. 

While we have emphasized the results here as predictions for the gauge theory correlation 
functions in an excited state at strong coupling, we believe that some of these predictions can 
be tested. For instance, one can consider the 1/2 BPS states in M = 4 SYM. The geometries 
dual to these states are explicitly constructed in [17]. Within this class of geometries one 
can focus on those which are asymptotically AdS and calculate the bulk prediction for the 
light-cone singularities using the nature of bulk geodesies. To compare with the field theory 
result we need to know the behavior of four point functions, which can be analyzed in large 
iV perturbation theory (at weak coupling) by taking the probe operators to also be chiral 
primary. We hope to report on this in the near future. 

12 This is because of the state being dual to a configuration of finite ADM mass M and the fact that the 
bulk Newton's constant scales like iV 2 . 
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4 Eternal black holes in AdS 



In the previous section we discussed the behaviour of geodesies in spacetimes which are de- 
formations of AdS, focussing on globally static geometries. We will now consider spacetimes 
such as the eternal Schwarzschild-AdS black hole which is causally nontrivial and not glob- 
ally static. In this case we know that a large Schwarzschild-AdS black hole corresponds to a 
thermal state in the field theory. The analysis of geodesies will lead to interesting signatures 
of the black hole geometry which will be visible in thermal correlation functions. 



4.1 Schwarzschild-AdS black hole: Structure of spacelike and null geodesies 

We will first describe the structure of null and spacelike geodesies in an AdS black hole 
background. Again focusing attention on 5-dimentional AdS, the black hole metric is given 
by 

ds 2 = -f(r)dt 2 + -^-dr 2 + r 2 dn\ (4.1) 
fir) 

with 

f(r) = r 2 + 1 - £ . (4.2) 

The mass of the black hole is proportional to /i, and correspondingly the horizon radius r + 
and the inverse Hawking temperature (5 are given by 13 



VI + 4/x - 1 Q _ 27rr H 



2 = v p = + / 4 _ 3 ) 

+ 2 ' 1 2r 2 + + l V 1 

The equations for general geodesies are given as: 

i = ~7~T\ i <P=^, and r 2 = a 2 - V eff (r) , (4.4) 
j{r) r z 

where 

K//(r) = ~/(r)+^ (4.5) 

and n = 1,0, or —1 for spacelike, null, or timelike geodesies, respectively, and the dot denotes 
differentiation with respect to the affine parameter along the geodesic. 

Null geodesies in the black hole background (k = 0) are parameterized by the ratio 
a = j. Once again we will be interested to know which points are connected by bulk null 
geodesies and we will quantify this in terms of (At(a), A(p(a)). The radial equation is easy 
to interpret as a classical particle with energy a 2 moving in a potential V e ff(r) = (see 
Fig. 6). V e ff{r) has a maximum at r 2 m = 2/i with a value V e ff{r m ) = 1 + ^ and as r — > oo, 
V e ff — > 1. The null geodesies connecting boundary points are: 



13 Note that to describe a big black hole we need (X > |. Further, for the temperature of the system to be 
high enough so that the big black hole dominates the thermal ensemble, we require /i > 2, which corresponds 
to T = > For illustration, we often use an example with /.t = 1 below. 
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Fig. 6: The plot for V eff (r) = f(r)/r 2 for /z e {1, 1.5,2,- •• ,7.5,8}. 



1. Geodesies with a = 1, i.e., E = J, stay at constant r = oo; these are null geodesies of 
the boundary manifold (Type A). 

2. The null geodesies which pass into the bulk and come back to the boundary (Type B) 
exist only for a G (1, a ), where 

a 2 = V eff (r m ) = l + j- (4.6) 

As a —» 1, we find that (At, Aip) — > (it, it). It is also clear from Fig. 6 that as a — > ao, 
At (a) and A(p(a) should go to infinity since it takes infinite affine parameter time to 
reach the turning point. 14 One can further check that 

dt dip dip , . , . „ s 

a— = -f = (4-7) 

da da dt 

From this we conclude that as a — > ao, 

At{a) w — Acp(a) -> oo (4.8) 

Note that since ao > 1, the above equation implies At < Aip. This remains true 
for general values of a G (1, ao) as is clear from Fig. 7, which is a parametric plot of 
Aip(a) — At(a) by varying a. Note that since <p is periodic, there is no violation of 
causality here. 



1 Or said differently, as a — > ao, the geodesic can go around the circular orbit at r m = many times 
before re-escaping to the boundary. 
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Fig. 7: End points of null geodesies in a black hole background. The black line is At = Atp, which is 
the set of end points of boundary null geodesies. The blue and red curves give the end points 
of bulk null geodesies for a black hole with \i = 2 and fi = 1, respectively. It is obtained by 
varying a G (1, cto). As a — > 1, (At, Atp) — > (vr,7r) and as a — ► arj, At and A<^ approach 
infinity with their ratio given by 1/ao- Note that Ap(a) > At (a) for a G (l,ao). 

In order to better understand the set of boundary points that can be connected through 
the bulk, we should also look at spacelike geodesies, which are discussed in detail in Ap- 
pendix C. 

4.2 New singularities for thermal YM theories on S 3 

We now consider gauge theory implications of the null geodesies in the black hole background. 
As discussed in Section 2, we will assume that null geodesies lead to the singularities of 
boundary correlation functions. 

At a spatial point Atp = <po, boundary null geodesies will give rise to singularities at 



These are singularities of vacuum correlation functions and we expect they survive at finite 
temperature as well. Type B null geodesies which go into the bulk will give rise to new 
singularities at 

At! = v?o + 2tt h(tp ), At 2 = At x + 2tt b 2 (tp ), ■•• , At n = At n _i + 2vr b n (tp ), • • • , (4.10) 

where At n , n = 1, 2, • • • are obtained by the intersection of the red line with vertical lines 
Aip = 27m + tpo in Fig. 7. Here b n (po) is some number lying between 1/ao and 1 and is 
smaller than 1. As n becomes large, 



At = <po, po + 2vr, p + 4tt, ■ • ■ 



(4.9) 
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independent of the value of (fo. The locations of these singularities are temperature-dependent, 
since «o an d b n depend on the temperature of the system (through //). Since all b n are smaller 
than 1, the singularities in (4.10) are generically distinct 15 from those in (4.9). Note as we 
increase the temperature of the system (i.e., increase the black hole mass /it), ckq — > 1 and 
the two sets of singularities merge in the fi — > oo limit. In the high temperature limit, S 3 
effectively decompactifies to R 3 , and as we discussed earlier there are no bulk null geodesies 
in this limit (see Appendix D) when we insist on maintaining the Poincare symmetry. 

Since "new" singularities arise from null geodesies which go into the bulk, they encode 
information regarding the bulk geometry. That singularities exist for At — > oo is a reflection 
of the existence of a (unstable) circular orbit at r m = a/2 /i in the bulk geometry. The time 
difference At n — At n _i between nearby singularities for n large is simply given by the time 
that it takes a null geodesic to transverse the circular orbit (which is 27r/ao). More generally, 
denoting the locations of the singularities by the curve At(Aip), from the second equation 
of (4.7) it follows that the slope gives 1/a = J/E of the corresponding bulk geodesic. 

Also note that in cases where the bulk geometry deviates from that of the Schwarzschild- 
AdS geometry in the vicinity of the horizon, such as the smooth microstate geometries 
discussed in the context of D1-D5 systems (cf., [18,19] for reviews and references), the gauge 
theory should in principle be able to detect even slight deviations which cause a slight shift 
in the radius of the null circular orbit. This is facilitated by the fact that the slope can 
be measured arbitrarily precisely for large enough separations At and Atp. 

Our conclusion above was obtained in the supergravity limit, i.e., large N and large 
A limit. It would be interesting to understand whether the "new singularities" at finite 
temperature arise when one departs away from the limit. In particular, it would be interesting 
if one could find independent arguments for their existence in the gauge theory. 

5 Horizon formation from gauge theories 

In the preceding sections we have focused on static 16 bulk geometries and indicated how to 
extract certain details about the geometry from the boundary correlators. In this section we 
finally apply the relation between null geodesies and singularities of boundary correlators to 
the fully dynamical case of gravitational collapse. Previous studies of black hole collapse in 
the AdS/CFT context are [20-23]. 

The basic question we want to ask is: Can we see horizon formation directly in the gauge 
theory? Since in the bulk, horizon formation is a sharply- localized event (despite the fact 

15 At some special values of <p and t, it is possible for them to coincide. 

16 Although the eternal black hole spacetime is not globally static, the method sketched above only probed 
the static part of the geometry. We comment on extracting the details of the dynamic spacetime inside the 
horizon in Appendix C. 
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that we need to know the full future evolution of the spacetime to actually find the horizon), 
we would expect that there will correspondingly be some sharp feature in the gauge theory. 
We will start by considering the simplest toy model, which is that of a collapsing spherical 
null shell. However, our main result should be applicable to general gravitational collapse, 
as will be clear below. 



5.1 Null geodesies in a gravitational collapse 

Consider a null spherical shell in d- dimensional asymptotically AdS spacetime. We can take 
the metric inside and in the past of the shell to be pure AdS, and the metric outside and to 
the future of the shell to be Schwarzschild-AdS. The metrics for the interior and the exterior 
of the shell can be written as 

dv 2 

ds 2 = -f in ,out(r) dt 2 in out + — + r 2 dVL 2 d _ 2 (5.1) 

Jin,out y ) 

where 

fin{r) = r 2 + 1 (5.2) 
fout{r) = r 2 + 1 - ^3 • (5.3) 

The mass of the black hole is proportional to fi. As in the previous section, we will use 
r + to denote the horizon radius and the surface gravity k of the black hole is given by 



27T 1 df out 



(3 2 dr 



(5.4) 



r+ 



with (3 the inverse Hawking temperature. Note that the t coordinate jumps across the shell 
with the jump approaching zero at the boundary. In contrast, the r coordinate is physical 
(since it measures the proper area of the spheres in a spherically symmetric spacetime) and 
therefore varies continuously across the shell. 

The Penrose diagram for the collapse is sketched in Fig. 8. There is only one asymptotic 
region; below the shell, the spacetime is pure AdS, whereas above, it is exactly Schwarzschild- 
AdS. 17 The two important times, which are labeled in Fig. 8 (a), are the time of the creation 
of the shell, t s , and the time of the horizon formation. In global coordinates, the latter 
occurs in the bulk at the origin r = at t = tn] however, this need not correspond directly 
to the CFT time, since we're describing an event at r = rather than at r = oo. Instead, 
the important CFT time related to the horizon formation is given by t = th (labeled on the 



1 The singularity is drawn with a tilt for more realistic representation of the Penrose diagram in d > 3 
dimensions (the shape of the singularity depends on d and r + ); however, here the "Penrose diagram" is to 
be treated as a sketch for ease of visualisation rather than an exact causal diagram. 
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(a) (b) 

Fig. 8: Penrose diagram for a collapsed black hole, (a) shows the times for the shell creation and 
horizon formation; (b) in addition illustrates three distinct regimes of radial null geodesies. 
(Note that the last (upper right) part of geodesic III is only relevant for nearly-null spacelike 
geodesies which can bounce off the singularity. Also, to avoid cluttering the diagram, the 
starting time U and ending time t are labeled only for geodesic I.) 
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boundary), which is the time from which a radial null geodesic would have to start in order 
to hit the horizon formation event (r = 0, t — tu)- 

Since this is a time-dependent geometry, the behavior of a null geodesic depends on the 
time ti of the initial point. Radial null geodesies therefore separate into three qualitatively 
distinct classes, indicated in Fig. 8 (b): 

I. For ti <t s — ir, the geodesic is given by that of the pure AdS, and t < t s . 

II. For t s — 7r < ti < th, the geodesic starts out in AdS, crosses the shell at some point 
outside the horizon, and continues in Schwarzschild-AdS, reemerging to the boundary 
at t Q G (t s , oo). 

III. For th < ti < t s , the geodesic starts out in AdS but crosses the shell at some point 
inside the horizon, and therefore hits the singularity (the last part of the geodesic 
III sketched in Fig. 8 (b) pertains to nearly null spacelike geodesies which can bounce 
off the singularity). For ti > t s , the radial null geodesic starts out in Schwarzschild- 
AdS and therefore crashes into the singularity; whereas for the nearly-null spacelike 
geodesies, this case is already covered above by reversing the orientation. 

We discussed geodesies in AdS and Schwarzschild-AdS spacetimes in previous sections; 
here we match them together to analyse the behaviour in the collapse geometry. We focus 
on null geodesies of Type B, with ti G (t s — ir,th); this will be the interesting regime where 
we can probe the event horizon formation. Note that th G (t s — ir,t s ), with th occurring 
earlier for larger black hole. In particular, using ray tracing within the AdS spacetime one 
can show that the horizon formation time is related to the shell creation time t s and its mass, 
or equivalently r+, by (see Appendix E for a derivation) 

t h = t s - 2tan~ 1 r+ . (5.5) 

Let us first consider radial null geodesies. From the Penrose diagram, it is clear that as 
ti — > th from below, t — > oo, while for ti > th, the null geodesic falls into the singularity and 
does not come back to the boundary. Since this will constitute the sharp CFT signature of 
the horizon formation event in the bulk, let us examine this feature in greater detail. The 
manner in which t Q diverges as St = tj L — ti — > + can be obtained by following Hawking's 
ray tracing argument [24], which gives 

t « -- log <ft + constant <ft ~ e~ Kto (5.6) 

K 

where k is the surface gravity of the black hole. 18 To see (5.6), consider a point P + on the 
horizon just outside the shell. Let v^ 1 to be the vector tangent to the horizon at P + and u 11 

18 Intuitively Eq.(5.6) reflects the fact that the redshift experienced by the outgoing null ray increases 
exponentially with t Q as e Kt ° when t — » oo. 
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be the null vector directed normally outward from the horizon normalized so that u ■ v — 1. 
The vector at P + which connects the horizon with the null geodesic can thus be written as 
ew' 1 for some e. As t — > oo, e can be expressed in terms of t a as 

e^Ae~ Kt % t a -> oo (5.7) 

with A some constant independent of t . Now parallel transport and across the 
shell to the point P_ on the horizon just inside the shell. 19 e remains invariant under this 
transformation. At P_ one can express e in terms of St as 

e^BSt (5.8) 

with some constant B. Comparing (5.7) and (5.8), we obtain the desired relation (5.6). 

Note that the leading behavior (i.e., the logarithmic dependence) in (5.6) is rather ro- 
bust, independent of the collapsing configuration and the type of black hole formed (charged, 
rotating). In Appendix E we also give alternative derivation of (5.6) for some explicit exam- 
ples. While our arguments so far were of a qualitative nature, we can make rather precise 
predictions for where the singularity in the collapse state correlators should appear. This 
is carried out in some detail in Appendix E, where we calculate t Q in terms of U and the 
parameters describing the state, t s , r + , and k, to be 

t = t s + — ( (1 + r + ) J - (1 + r\) tan" 1 — ^ + r + tanh^ 1 — ) (5.9) 
r + k \ 2 1+7+ r + J 

where r c denotes the crossing radius where the geodesic intersects the shell, r c = tan 

To summarize, we find that radial null geodesies are sensitive probes of horizon formation. 
The infinite redshift at the horizon translates into a sharp time scale t^ (5.5) at the boundary. 

Now let us consider non-radial null geodesies. A natural question is whether there exists 
an analogous value th{a) for each 20 a = E/J so that as U — > th(cx), the corresponding 
t — > oo. The answer turns out to be yes for a range of values of a, but for a different 
reason from that for the radial geodesies. By the very nature of an event horizon, we expect 
that the only null geodesic which truly samples the event horizon formation is the radial 
one; any other null geodesic which approaches the horizon non-radially cannot escape to 
the boundary. In Appendix F, we describe how to compute th(a) for general a. We show 
that all geodesies with angular momentum for which t D diverges take that long to reach the 
boundary not because of being trapped near the event horizon, but rather because of circling 
around a null circular orbit at a given finite distance outside the horizon. 

It is easy to see that th(a) depends non-trivially on a; the presence of angular momentum 
causes the geodesic to sample a different region of the spacetime geometry. In the limit 

19 Note that m m and are continuous when crossing the shell. 

20 Since the shell geometry is not static, E jumps across the shell; in defining a we use the initial E i.e., 
its value in pure AdS. 
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a — > 1 (maximal angular momentum) we expect that the geodesic stays arbitrarily close to 
the boundary, so that it doesn't sample the black hole geometry, and in particular there can 
be no divergence in t a for any finite ij. In the opposite limit of zero angular momentum, 
we have already seen that t h (oo) is given by the expression for t h for the radial geodesic 
(5.5). The result of Appendix F is summarized in Fig. 9, which plots th{ot) as a function of 
J ~ 1/a. In particular, we note that: 

th(a) — > th = t s — 2 tan _1 (r i ) , a — > oo 

(5.10) 

t h (a) ->■ t s , a — > a 

while for 1 < a < a^, t never diverges, where ao was introduced in (4.6). Recall that the 
Schwarzschild-AdS geometry contains a null circular orbit with a = a at r = y/2 Ji. As 
Fig. 9 demonstrates, th{a) is not a monotonic function of 1/a. In particular, the minimal 
value of th(a), which we will denote as t c , is smaller than th, the horizon formation time. t c 
is the boundary time scale at which the null circular orbit of the newly formed black hole 
geometry is first probed. 

t h (a) 
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Fig. 9: Variation of th(a = j) as a function of angular momentum J in the thin shell spacetime (5.1) 
(solid line). The intersection of the curve with vertical axis gives th, while the minimum gives 
t c . Red dots are data extracted from Fig. 11, which were obtained by numerical integration of 
the geodesies in the Vaidya spacetime (5.12), (5.13) with v s = 0.001. 



In Fig. 10, we plot the endpoints of various geodesies on the At — A(f plane 21 for a fixed 
initial time ij G (t s — vr,t s ). Note that in contrast to the analogous plots in earlier sections 
for other geometries, the background here does not have a time translational symmetry. For 
ease of visualization we unwrap the tp direction and at the same time compactify both the 

21 Since a is the only continuous free parameter describing the null geodesies, this will produce a curve of 
endpoints on the At — Aip plane. 



23 



At and Aip so that we can examine the full At — Aip plane. For each fixed ti, we plot the 
end-point curve, color-coded by a. For ti < t s — ir, the endpoints would all clump into the 
single point (Aip, At) = (vr,7r). As we increase ti so that the geodesies start sampling the 
shell, the endpoints begin to spread in the manner shown on the left in Fig. 10, into a cusp 
similar to the star geometry (cf. Fig. 5). As ti reaches the minimum of the curve th(ot) shown 
in Fig. 9 (i.e., t — > t c ), the cusp extents to (Aip, At) = (oo, oo) (upper right corner of the 
plots) with 

At(a) » — A<p(a) -> oo . (5.11) 

Increasing ti further, the cusp still reaches (Aip, At) = (oo, oo) with (5.11), but now for two 
different values of a, given by solution to th(a) = ti. At the same time, At for the radial 
geodesic increases, diverging as ti — > t h . This will appear in Fig. 10 as the left end point of 
the upper branch moving to (it, oo). The right plot in Fig. 10 is taken for ti slightly smaller 
than th- Note that equation (5.11) follows from the fact that as ti —>■ th(a), the geodesic 
goes around the circular null orbit with a period — infinitely many times. 
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Fig. 10: Plot of (Aip, At) for 3 values of U in the thin shell geometry (same configuration as in Fig. 9 
and Fig. 11(a)), for various values of a G (1, oo). The horizontal and vertical axes show Ap 
and At respectively, rescaled by tan" 1 to fit the full plane. The left endpoint of the upper 
branch corresponds to a = oo (radial geodesic), while the left endpoint of the lower branch 
corresponds to a = 1. (a) U = —2.8 < t c « —2.72; (b) ti = —2.73 only slightly smaller 
than t c ; (c) U = —2.67 slightly smaller than th ~ —2.66. Here t c is the minimal value of 
the curve in Fig. 9 and th is intersection of the curve in Fig. 9 with the vertical axis. As 
ti — ► t c the cusp should reach (oo,oo) (upper right corner) and as ti — > th, the left endpoint 
of the upper branch should reach (n, oo) (the upper left corner). Note that as apparent from 
Fig. 11, the endpoints vary rather sharply with ti for ti th(ot), so it is difficult numerically 
to sample the large values of Aip, At. 

One can generalize the sharp shell to a smeared-out version which is physically more 
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realistic. To that end, we consider a Vaidya spacetime of the form 

ds 2 = -f(r,v)dv 2 + 2dvdr + r 2 dflj (5.12) 

with f(r,v) smoothly transitioning between (5.2) at v — > — oo and (5.3) as v — ► +oo; a 
convenient form to use is 

u fl +tanh^\ 

/(r, v) = r 2 + 1 - £ f (5.13) 

In particular, the shell is inserted at v = = t s + tt/2, and has 'thickness' v s . As v s — > 0, we 
recover the collapse spacetime (5.1) written in ingoing coordinates. 




Fig. 11: Plot of At as a function of tj for non-radial null geodesies in the Vaidya metric with the 
effective "thickness" of the shell given by (a) v s = 0.001 and (b) v s = 1. The various 
curves differenct angular momenta, J = 0.01, 0.1, 0.2, 0.3, .... 



Fig. 1 1 shows the value of At = t Q — tj for a set of non-radial geodesies in the spacetime 
(5.12) with various values of a, all starting at the initial time tf, in (a) the shell is very thin 
whereas in (b) the shell is considerably smeared out. The value of th(a) can be read off from 
the plot; it corresponds to the value of at which At diverges for a geodesic with that a. 
For the thick shell in Fig. 11 (b), we see that these times are well-separated, and increase 
monotonically with increasing J/E (decreasing a). 

This illustrates that the sharp signature of the horizon formation event is not limited to 
the thin-shell collapse (bulk spacetimes with a null shell sharply separating pure AdS from 
Schwarzschild-AdS): For any spherically symmetric spacetime in which an event horizon 
forms at a time tn at the origin, reachable by an ingoing radial null geodesic starting at 
the boundary at t^, we can find th by exactly the same method. In particular, our method 
constitutes following the red (leftmost) curve in Fig. 11 (b). 

To summarize, whereas only the radial null geodesic starting at th samples the horizon 
formation event, the non-radial null geodesies provide further details about the geometry, in 
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particular the null circular orbit of the newly formed black hole geometry. Furthermore the 
variation of the curves At(A<p) with respect to t$ provides spatio-temporal information both 
dynamically around the horizon formation time and spatially in the vicinity of the horizon. 

5.2 Gauge theory signatures 

We can now translate the behaviour of null geodesies in the bulk collapse spacetime to 
predictions regarding boundary correlation functions. Our discussion below should apply to 
generic highly excited states in the boundary theory 22 . 

Consider the CFT correlators of the form (0(ti,Q) 0(t , —Q)) shcU evaluated explicitly 
in the state of the shell. Modelling the shell as being created by an operator 5 inserted in 
the boundary at t = t s , we have 

( o{t , n') o(u, n) ) shcll = (o| o{t a , n') o(u, n) |o> , u < t < t s 

= (0| 0(t , 0') Ofo, Q) 5(t s ) |0) , U<t s <t (5.14) 

= (o| s\t s )o{t , n') o(u, n) s(t s ) |o) , t a <u<t 

By our previous arguments, we would expect this correlator to become singular when the 
insertion points are connected by a null geodesic. More explicitly, 

I. When t < t s , the correlator is simply the standard vacuum two point function. The 
singularities are given by the usual light cone singularities of the boundary theory. As 
is clear from geodesic I in Fig. 8(b)), all null geodesies lie entirely in the AdS region. 

II. When ti < t s < t , the correlator is given by the second line of (5.14). From our 
discussion of null geodesies in last subsection, the pattern of bulk-cone singularities for 
this correlator has a rich structure, which reflects the formation of the event horizon 
and the null circular orbit in the bulk. 

Suppose experimentalists in the boundary theory are able to measure the observable 
shell f° r an ^ an d to satisfying ti < t s < t Q . Then by carefully 
plotting the locations of singularities of the correlator in the At — A(p plane for a given 
ti, they should recover various plots of Fig. 10. In particular two distinct time scales 
emerge by comparing the pattern of singularities for different U. The first is the time 
t c when the cusp in Fig. 10 reaches (oo, oo). This is the time scale that non-radial null 
geodesies originating from the pure AdS geometry start being trapped by the circular 
orbit of the newly formed black hole. The slope of the line reaching the infinity gives 
the period of the orbit. The second is the time t^ that the left end of the upper branch 

22 By highly excited states we mean states of energies of order cN 2 , with N- independent constant c suffi- 
ciently big. By generic states we mean states which are generic superpositions of energy eigenstates. Note 
that energy eigenstates do not give rise to the desired time dependence. 
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(which corresponds to radial geodesies) reaches (it, oo). This is the time of horizon 
formation as probed by a radial null geodesic. The bulk time of horizon formation in 
global AdS coordinates is given by tjj = th + |. 

III. When t s < ti < t Q , the correlator is given by the third line of (5.14). In this case, all 
singularities of the correlator are the same 23 as those of a finite temperature correlation 
function with the temperature T given by 

1 9\o g Q(E) 

T~ dE {bAb) 

where E is the energy of the shell state and Q(E) is the density of states of the CFT. 
This follows from that all type E geodesies with t s < t { < t are identical to those of 
an AdS black hole. 



In the above we discussed the pattern of bulk-cone singularities of (5.14). There might 
also be more subtle signals, on a secondary sheet (so we would find the singularity in the 
correlator only after a suitable analytic continuation) when the insertions are separated by a 
spacelike geodesic which is arbitrarily close to being null as illustrated for the eternal black 
hole singularity in [6]. In other words, the corresponding geodesic passes through the black 
hole and bounces off the singularity 24 , see e.g. geodesic III of Fig. 8(b). Such a situation 
could arise when ti > th (which, if we define ti < t Q , also implies that ti < t s < t Q ). We will 
discuss it in more detail in Section 5.3. 



5.3 Signature of the black hole singularity? 

In addition to the null geodesies considered in Fig. 8, it is also of interest to consider some 
spacelike geodesies. In particular, spacelike geodesies have the important feature that they 
'bounce off' the black hole singularity [6], and thus could provide valuable information about 
the black hole singularity. There are two other types spacelike geodesies in the collapse 
background which are of particular interest (see Fig. 12): 

1. ti — > th from above and t D — > oo. As can be seen from Fig. 12a, the geodesic passes inside 
the horizon and bounces off the singularity. We will calculate below the dependence of 
t Q on 5t = th — ti. 

2. As ti — > t s , t Q — * t s . The corresponding geodesic is sketched in Fig.l2b. One expects 
that t — t s ~ (t s — ti) 1 , where 7 should depend on how strongly geodesies are repelled 
from the singularity. 



23 The correlator itself does not necessarily coincide with a thermal correlation function. Only the singu- 
larity structure does. 

24 Howcver, note that there is no longer a single radial null geodesic connecting the insertion points. 
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(a) 



(b) 



(c) 



Fig. 12: Penrose diagram for a collapsed black hole 



Note that spacelike geodesies are parameterized by E and J, but here we set J = to focus 
on radial geodesies, and moreover we take the infinite E limit, so that their trajectories 
converge to the radial null geodesies, except for the important point that they 'bounce off' 
the black hole singularity As before, we can find geodesies in the full collapsing geometry 
by patching the geodesic segments in the respective regions; here in addition to AdS (region 
1 in Fig. 12c) and the Schwarzschild-AdS region outside the horizon (region 2 in Fig. 12c), 
we now have to consider Schwarzschild-AdS inside the horizon as well (region 3 in Fig. 12c). 
This is presented in detail in Appendix E; here we only quote the main results. 

1. As ti — > t^, t Q — > oo in the same manner as for — > t7 (cf. Eq.(5.6)), namely: 

St = U - t h ~ e~ Kto (5.16) 

2. Let t s — ti = e. Then expanding t Q — t s to third order in e, we find that 

t -t s ~e 3 . (5.17) 
with a coefficient which depends on the parameters of the final black hole - see (E.16). 

As argued around (5.6), it is clear that the behaviour (5.16) is robust and is expected 
to hold in any spacetime dimension. The behaviour of the geodesic which bouces off the 
singularity in Fig. 12b, and in particular the exponent 7, depends crucially on the dimension; 
7 = d — 2 in AdSd. Intuitively this is due to the fact that the curvature of the singularity 
increases with spacetime dimension. 
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6 Discussion 



We have argued that the AdS / CFT correspondence implies the existence of "bulk-cone sin- 
gularities" in boundary theory correlation functions, which lie inside the light-cone. The 
pattern of singularities can in turn be exploited to extract information about the bulk geom- 
etry. Our argument revolved around the essential point that Green's functions for quantum 
fields in curved space are singular when the arguments are null separated. While the geodesic 
approximation was necessary to derive the conclusion, we should point out that one expects 
the result to be true as long as it makes sense to talk about the spacetime geometry. In 
other words, we expect validity in a 1 expansion, but not necessarily in g s expansion (where 
it is plausible to have a breakdown of the geometric picture). 

We have demonstrated the pattern of singularities in the field theory correlators for 
several distinct scenarios: in radiation star, eternal Schwarzschild-AdS black hole, and the 
geometry of a collapsing shell. We have seen that small deformations of bulk geometries 
can in principle produce significant differences in the structure of singularities which are 
discernible on the boundary. In this context it is also satisfying that the vacuum state of 
the CFT, corresponding to pure AdS spacetime, has a distinguished role, as evidenced by 
the theorem of [3]. Further, we are guaranteed that bulk-cone singularities for correlators 
evaluated in excited states occur at later times compared to the vacuum correlators, as 
necessitated by causality. 

The novel application of our proposal was to the case of dynamical spacetimes. We con- 
centrated on the simple scenario of a null shell collapse in AdS/CFT. One of the fascinating 
aspects of the singularity structure of the correlations is that they encode the spacetime 
event corresponding to event horizon formation in the bulk unambiguously. The charac- 
teristic signature we discussed is quite robust and can be easily distinguished from other 
singularities in the correlation functions. Further, the detailed pattern of singularities car- 
ries information about the geometry of the black hole that is formed. By using the non-local 
nature of the correlation functions we are able to extract information about the details of the 
bulk spacetime, while, in contrast to the recent discussions of probing the black singularity 
in AdS/CFT [6,7], we are not reliant on analytic properties of the CFT correlators. 

It would be very desirable to test our predictions of the existence of the bulk-cone singu- 
larities directly within the boundary theory. Note that the classical limit corresponds to the 
large N and large 't Hooft coupling limits of the boundary theory. It would be interesting 
to understand whether the appearance of bulk-cone singularities for correlation functions in 
excited states is an artifact of the large N and/or large 't Hooft coupling limits, or a generic 
feature for field theories on compact spaces. 25 Further, we expect that the semiclassical pic- 
ture that we used in the bulk to be valid only in the large N limit - at any finite N we will 

25 Note that bulk-cone singularities do not arise for non-compact boundary as shown in Appendix D for 
states preserving the boundary isometries. 
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have contributions from other semi-classical geometries that are sub-dominant saddle points 
of the quantum gravity path integral. In these circumstances it is not clear whether there is 
a precise meaning to the bulk light-cone. Nevertheless, we should note that the bulk-cone 
singularities discussed here are eminently trackable in the field theory even beyond the planar 
limit of large N. Understanding the role of these in some geometric terms is an interesting 
challenge. 

An intriguing avenue to explore is to extend these considerations to "microstate geome- 
tries" proposed for black hole spacetimes. In situations where we have an explicit map 
between the spacetime geometries and the field theory states such as two charge D1-D5 sys- 
tem [18,19], the 1/2 BPS sector of M = 4 SYM [17], one can ask whether the pattern of 
singularities differs significantly between different states. The general arguments regarding 
the correlation functions of typical states that comprise the black hole entropy in [25-28] 
would seem to suggest that the structure of the singularities is not capable of discerning the 
fine distinctions between the geometries. This is however counter to the gravitational intu- 
ition: naively one expects geometries that differ from each other to have differing behaviour 
of Green's functions, especially with regard to the location of the singularities. A potential 
resolution of this discrepancy is that the geometric picture is incorrect for the "microstate 
geometries" as one has to consider a suitable wavefunction over the quantum moduli space 
of solutions. This is an interesting question that deserves to be explored further. 

Another interesting generalization of our set-up would be to model the AdS/CFT analog 
of Choptuik scaling [29] . Given that we have a specific prediction for the field theory signature 
of black hole event horizon formation, one can ask whether the critical behaviour observed 
during the collapse has a field theoretic image and map out the details of the correspondence 
in this case. 
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Appendix A: Pure AdS 



For ease of comparison of our results on the star and other geometries in what follows, 
we first summarise the basic facts regarding geodesies in AdS, concentrating on null and 
spacelike geodesies. 



A.l Geodesies in AdS: summary 

Consider d + 1-dimensional AdS in global coordinates. The metric can be written as 

f{r) (A.l) 

f(r) = 1 + r 2 . 

The geodesic equations are as in (4.4), with f(r) being that given in (A.l). 

As we are primarily interested in the endpoints of null or spacelike geodesies (ultimately 
with diverging regularised proper length). To that end, the useful quantities are the tem- 
poral and angular separation of the geodesic endpoints, which we denote by At and Aip, 
respectively. Let us first focus on spacelike geodesies, since as we will see, these limit to null 
geodesies. We can integrate the geodesic equations to obtain: 

At = - + sin" 1 — (A.2) 



2 V e ~ L 

n 1 {E 2 -J 2 + 1\ 
Aip = - + sin" 1 — (A.3) 



2 ' V e " J 

e~ L = yj (E 2 — J 2 ) 2 + 2(E 2 + J 2 ) + 1 (A.4) 

where L is the geodesic length after regularising away the divergent piece log R 2 C (Rc is the 
radial cut-off). Finally, another useful expression is the minimum r value reached by the 
geodesic r min 

2r 2 min = -(E 2 -J 2 + l) + e- L (A.5) 
From the expression for e~ L it is clear that we have a large negative length when 

1. E — > oo with J fixed =^ L ~ — log-E 2 

2. J — > oo with E fixed => L ~ — log J 2 

3. J ~ E — ► oo L~— logE 
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Now consider two boundary points A and B are separated by boundary coordinate dis- 
tance At = 7r — e t and Ay? = ir — e v . We can easily see that for any spacelike separation 
< e^ip < 7r, there is a unique geodesic joining A and £>; specifically we can calculate E and 
J for the geodesic given e t v : 

£ = sinc « (A.6) 

cos — cos e 4 

J = Sin ^ (A.7) 
cos e v — cos e t 

and its length is 

L = — log 2 + logf^cose^ — cose t ) (A. 8) 



A. 2 Spacelike vs. null geodesies 



Null geodesies in AdS only connect points with At = Ay? = it, irrespective of the parameter 
E/J. While this is a well known fact for radial null geodesies (for which J = 0), it is true that 
even geodesies carrying angular momentum emerge at the anti-podal point on the sphere in 
exactly AdS time. As mentioned in Section 2 these null geodesies endpoints being coincident 
with the endpoints of purely boundary null geodesies, we do not see extra singularities in 
the correlation functions. 

A naive puzzle about the null limit arises when we consider the difference between the 
angular separation of the endpoints A<p for spacelike and null geodesies. Spacelike geodesies 
in AdS allow any < Ay? < 7r; indeed for E = (constant t slice) the geodesies plotted on 
the Poincare disk are circular arcs with diameter related to the angular momentum. On the 
other hand, null geodesies always reemmerge at the anti-podal point Ay? = it. How do the 
geodesies change their behaviour as we take E — > oo at fixed J? Fig. 13 illustrates this limit. 
For any fixed finite J, as E — > oo, we see that the geodesic indeed converges to the radial 
null geodesic and Ay? — > ir. 

Nevertheless, the full set of singularities of the boundary correlation function (2.10) is 
larger than those captured by the points with At = Ay? = ir. It is interesting to ask whether 
there are any bulk geodesies (with divergent regularized proper length) that connect points 
with At = Ay? ^ 7r. From (A. 4), one concludes that the regularised proper length along a 
spacelike geodesic diverges not only when E — > oo, but also when J — > oo. Does this mean 
that this is an equally good null limit, and that therefore Ay? = 7r? The answer is no: as 
we can see from (A. 3), as J — > oo at fixed E, Ay? — > 0. The geodesic has vanishing proper 
length because it is so short rather than because its tangent vector is null. On the boundary, 
this corresponds to the usual divergence of correlator of operators inserted at the same point. 

However, it turns out that there nevertheless are nontrivial almost-null spacelike geodesies 
with At = Ay? 7^ 7r. Consider spacelike geodesies in AdS parameterized by E and J with 
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Fig. 13: Spacelike geodesies in AdS, projected to a constant t slice (left) and the t — r plane (right), 
for six different values of energy (E = 0, 2, 5, 10, 30, 80). For each E, geodesies with different 
values of angular momentum ( J = 0, 1, 2, . . . , 10) are plotted. The bold circles represent the 
boundary at tanr = ^ The dashed vertical lines correspond to the origin r = and the 
bold vertical lines to the boundary at tanr = ^. The range oft plotted is (0,tt). 



the following constraint: 

J 2 = E 2 -aE (A.9) 
for some real number a. Upon taking the limit E — > oo, we obtain 

At -> g + sin' 1 ( - a ^ (A.10) 



2 



2 Vv^n, 

Clearly, these are equal, and take the full range < At = Ay? < n for — oo < a < oo. 
Further, we can check that this family of geodesies indeed converge onto the boundary null 
geodesies. Fig. 14 demonstrates these features. We see that as E increases, the geodesies 
accumulate at a finite value of At and A<^, in this case a being 1, | + sin -1 (^^j ~ f • 

To summarise this discussion, in pure AdS^+i spacetime, null geodesies connect points 
on the boundary which are anti-podally located on the S d_1 at a time separation At = ir. 
Spacelike bulk geodesies on the other hand connect points with At = Aip ^ ir as illustrated 
by (A. 10) and (A. 11), respectively (or more generally, (A. 2) and (A. 3)). 
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(a) 



(b) 




Fig. 14: Spacelike geodesies in AdS, projected to (a) a constant t slice and (b) the t — r plane, 



Appendix B: Radiation star in AdS 

In this Appendix, we construct the geometry corresponding to a gas of radiation in AdS 
and study its properties. As sketched out in Section 3, constructing the solution entails 
specifying the form of the stress-energy-momentum tensor and solving the Einstein's equa- 
tions. 26 After filling in the gaps in the presentation of Section 3, we extend our study of the 
star geometry by analysing spacelike geodesies, and then discuss the observable differences 
between the star geometry and that of a black hole. 

B.l Construction of the star- AdS spacetime 

As discussed in Section 3, one can construct a simple model for a static, spherically sym- 
metric, asymptotically AdS "star" geometry by solving Einstein's equations with a negative 
cosmological constant and matter given by a perfect fluid stress tensor corresponding to 
radiation. The metric coefficients are determined by solving the Einstein's equation 



where for convenience we will set 87rG5 = 1 and A = —6 to set AdS radius to unity. 

The symmetries constrain the metric to take the form (3.2). One can then infer the 
equations (3.3) and (3.5) from the tt component of the Einstein tensor. The rr component 

26 Analogous calculations in 4 dimensions were discussed in [13], with similar qualitative results as we 
obtain below. 



for increasing values of energy and J 2 = E 2 — a E with a = 1 (in the infinite E limit this 
would yield At = A(p « ^). In (a) the boundary is indicated as the bold circle. In (b) the 
boundary at tan r = | is in bold and the dashed vertical line is the origin. 
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of the Einstein tensor yields 



L 
f 




(B.2) 



which can be easily integrated to get f(r) in terms of m(r) and p(r). Finally, using the 
angular part of Einstein's equation, or equivalently the stress tensor conservation, we can 
write (for general equation of state) 

dP If 

£ = + ^ (B.3) 

This can then be used to derive the system of coupled first order ODEs for m(r) and p(r) 
given by (3.6) and (3.7), respectively. The equation for f(r), (B.2), can be simplified to (3.4) 
using the equations for m(r) and p(r). 



B.2 Geodesies in the star geometry 

In Section 3 we examined some aspects of the star geometry by focusing on null geodesies; 
here we extend this analysis to spacelike geodesies as well. While these do not lead to 
singularities of the boundary correlators, they nevertheless reveal interesting points. Since 
the star in AdS can be vied as a deformation of the pure AdS geometry, it is particularly 
interesting to contrast the behaviour of spacelike geodesies in the star geometry with that 
in AdS. In Fig. 15 we show spacelike geodesies projected onto the constant t and t — r plane 
(these plots are to be compared with Fig. 13, where same conventions apply). As for null 
geodesies, we see a focusing effect, created by the star in the center of AdS. 

To see how the endpoints of these various geodesies compare, let us consider the values At 
and Aip for the various geodesies. Fig. 16 summarizes this. One interesting feature to note 
is that since the curves for various energies intersect each other, it is no longer true that for 
any two spacelike-separated points, there is a unique spacelike geodesic connecting them. In 
particular, unlike the pure AdS case, there is an open set of endpoints which can be reached 
by more than one geodesic. Even more remarkably, there exists a set of endpoints (lower 
branch of the thick pink curve in Fig. 16) which are connected by both a null and a spacelike 
geodesic! Also, as a corollary, it is no longer true that the endpoints of null geodesies lie on 
the boundary of the set of endpoints of spacelike geodesies. Correspondingly, the enpoints 
of the null geodesies reflect this behaviour, as can be seen easily in Fig. 5. While Aip slowly 
increases with increasing po, the endpoints never approach a straight line, as would be the 
case for the black hole geometry. This in turn means that the endpoints are less sensitive 
indicator of the star geometry than in the black hole case. 
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Fig. 15: Spacelike geodesies in star with po = 10 in AdS, projected to a constant t slice (left) and 
the t — r plane (right), for six different values of energy (E = 0,2,5,10,30,80). For each 
E, geodesies with different values of angular momentum (J = 0,1,2,... ,10) are plotted. 
The bold circles represent the boundary at tanr = |. The dashed vertical lines correspond 
to the origin r = and the bold vertical lines to the boundary at tanr = \. The range of t 
plotted is (0,1.1 At ). 



B.3 Comparison between star and black hole 

Let us now ask how does the geometry vary with po- Clearly, as po — ¥ 0, the geometry 
becomes that of pure AdS; hence the star's utility in studying small perturbations on AdS. 
One might naively expect that in the opposite limit, as we make the internal density very 
large, the star should start behaving more and more like a black hole. Unfortunately, this is 
not the will be explained below. 

First of all, to compare the star's geometry with that of a Schwarzschild-AdS black hole, 
the most sensible map of parameters is to identify the black hole's mass with the star's total 
mass, so that the asymptotic geometry matches. In other words, for the black hole, we use 
the metric of the form in (3.2), but with 

r 2 M 

f( r ) = h (r) = — + 1 , where M = lim m(r) (B.4) 

Now, consider the variation of the star's density and mass profiles as we increase po- This is 
plotted in Fig. 17. We see that as the star's total "extent" does not seem to change much with 
Po - the star is always confined within a size of order the AdS radius. This is of course due to 
the confining potential of the AdS geometry. Now, if the total mass could be made arbitrarily 
large, then for M > 1, the star would be confined within its own Schwarzschild radius, and 
should be viewed as a black hole. However, we see that while the density can be increased, 
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16: Endpoints of spacelike and null geodesies in the geometry of star with po = 10 in AdS. 
Spacelike geodesies for six different values of energy (E = 0,2,5,10,30,80) and a null 
geodesic (thick pink curve) is plotted for varying values of J - varying J traces out curves 
in the At — Aip plane; the top of the curves corresponds to small values of J, while in the 
large J limit, At and Aip both approach their starting value, namely 0. 



p(r) m(r) 




17: Density profile, and corresponding mass profile, for a star with various values of po (which 
can be read off from values of p(r = 0)). 
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the mass is bounded from above, and becomes largely independent of po after some point. 
This is illustrated more clearly in Fig. 18, where (a) the total mass, and (b) the fraction of 




Fig. 18: Star's total mass (a), and the fraction of the star's mass confined to within this radius, 
plotted as a function of po- 



the star's mass confined to within its effective Schwarzschild radius is plotted as a function 
of po- The left plot shows clearly that the total mass is bounded 27 from above by a rather 
small value: M < 0.775R. Curiously enough, the mass is not a monotonic function of the 
internal density. 28 In fact, when examined in more detail, M(p ) appears to exhibit certain 
self-similarity; however we will not detour into this intriguing observation further. We also 
find that the effective Schwarzschild radius (i.e., what would be the Schwarzschild radius for 
a black hole with the total mass plotted in Fig. 18a) is bounded by r + < 0.716i?, so the star 
could at best be compared with a small Schwarzschild-AdS black hole. Most importantly, 
Fig. 18b shows that the mass contained within the effective Schwarzschild radius is a small 
fraction of the total mass (less than 11%), so that the star does not approach a black hole 
in any regime. In other words, we should not expect to obtain behaviour characteristic of a 
presence of an event horizon for any value of po- 

We have seen that the star made up of a gas of radiation can never develop an event 
horizon, so it cannot look like a black hole, no matter how high its internal density po- m 
Section 3 we have further argued for the absence of null circular orbits, and emphasized that 
the corresponding geodesic endpoints demonstrate a clear difference between a star and a 
black hole. To illustrate this further, in Fig. 19 we plot the null geodesies, projected to a 
constant t slice, for increasing values of po. (This is the same plot as in the left plot in Fig. 4, 

27 In fact, this result is analogous to the similar result in the more familiar 4-dimcnsional, asymptotically 
flat static spherically symmetric spacetime, where given a fixed size i?* of a star (with any equation of state), 
the maximum possible mass such a star can attain is A/ max = 4i?*/9. [30] 

28 This might correspond to onset of some sort of radial instability, though there docs not seem to be any 
obvious pathology. Correspondingly, we'll continue to include such high-po solutions in our considerations. 
For the stars in 4-dimensional, asymptotically flat static spherically symmetric spacetime, [31] argued that 
an instability sets in at the turnover point where the mass function is non-monotone. 
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Po = 100 



Po = 1000 



p 



= 10000 




Fig. 20: Plots of V e ff(r) for J = 5,3,0.1, which correspond to red, green and blue curves respectively. 

The horizon is at r + = .78 (for /i = 1). For J < r + (e.g. blue curve), geodesies will fall into 
the horizon and will not come back to the same boundary for any E. Similarly for E > E max . 
Emax is the maximal value of V e ff. 

which had po = 10, redone for po = 100, 1000, and 10000, as indicated.) We see that even at 
very large po, there is no null orbit. Even though the bending increases with po, it does so 
ever more slowly, rather analogously to the time delay for radial geodesies (cf. Fig. 3). 



Appendix C: Eternal black hole in AdS 

We now discuss some interesting properties of spacelike geodesies in the Schwarzschild- 
AdS background. The geodesic equations are given in (4.4). The quantity of interest for 
spacelike geodesies is (At(E, J), A<p(E, J)) which depends on both E and J. 

The effective potential V e ff(r) given in (4.5) for spacelike geodesies has two real zeros 
at r = J and r = r + , respectively. For J < r+, V e ff(r) < for r > r+, so for any 



39 




Fig. 21: The end points of various spacelike geodesies in AdS BH (with \i = 1 and r + = 0.78). Black 
and red line correspond to boundary and bulk null geodesies respectively as in Fig. 7. (a): 
Various curves show how end points change as E is varied from E = to E max (J) with 
fixed J corresponding to J = 0.8, 1.2, 2.2, 6.5. The curve with larger J is farther away from 
the the lines corresponding to the null geodesies. Note that at E = 0, all lines start at real 
axis and as E — > E max (q), At and Aip go to infinity, (b): Various lines describe how end 
points change with J for a fixed a = E/J = 1.1, 1.08, 1.02, 1, 0.2. The line with bigger 
value of a is closer to the lines of null geodesies. J varies from some finite J c to oo. It is 
clear that as J — > oo all curves approach their corresponding points on the null curve. Note 
that a curve with < a < 1 will end up at the origin as J — > oo. A curve with a = 1 could 
end up anywhere in the diagonal line between the origin and (vr,7r) depending the value of 
E — J which should be fixed in the limit E, J — > oo. 



E, the spacelike geodesic will fall into the horizon 29 and will not come back to the same 
boundary. We are not interested in these geodesies. For J > r + , V e ff(r) has a maximum 
V C (J) = V e ff(r m ) at a value r m > r + (see Fig. 20). For E 2 > E^ ax (J) = V C (J), the geodesic 
will again fall into the horizon. Thus for a spacelike geodesic to come back to the same 
boundary, J and E should lie in the range J G (r + , oo) and E G (0, E max (J)) for a given J. 

For a fixed J > r+, at E = 0, (Acp(0, J), At(0, J)) = (A(p c (J), 0). A(p c (J) is a monotonic 
function of J. As J — > oo, Aip c ( J) — > j — > 0, while as J — >• r+, A(p c (J) — > oo logarithmically. 
As E varies from to E max (J), At(E, J) and A(p(E, J) appear to monotonically increase to 
infinity (see Fig. 21). As E — > E max (J), it can be shown that 

At(E, J) w -i- Ay?(£, J) - oo (C.l) 
s(J) 

where s(J) is given by 



•w = = + (c ' 2) 



2!) 



Subsequently it will either reach the other boundary or fall to the singularity. 
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Recall that r m (J) is the maximum of Vj(r). Note that r m (J = r + ) = r + and as J — > oo, 
v m {J) — ► (2/Lt)5. One can also check numerically that r m (J) is a monotonic function of J. 
From the second equality of (C.2), s(J) is then a monotonic function of J decreasing from 
s(J = r + ) = yJ2 + Jj- to s(J = oo) = a , where a was introduced in (4.6). 

From the above discussion we conclude that the end points of spacelike geodesies cover 
the whole region bounded by the red curve corresponding to the bulk null geodesies, the 
straight line from the origin to (it, it), and the horizonal axis in Fig. 7. It is important to 
emphasize that (f is a periodic variable, so in various figures, the fundamental region is the 
strip if G (— 7r, 7r]. As a result we conclude that 

1. For any end point (Atp, At) of a null geodesic (including boundary ones), there are an 
infinite number of spacelike geodesies ending at that point due to the periodicity of <p. 

2. The bulk and boundary null geodesies can also share their end points for special values 
of At and Aip. 

Finally let us comment on how to obtain null geodesies from spacelike geodesies by taking 
a limit. It is easy to see that we see that the equations for spacelike geodesies approach those 
of null geodesies in the limit of E, J — ► oo while keeping the ratio a = E/J fixed. More 
explicitly, the end point of a spacelike geodesic with < a < 1 will approach the origin 
in the J — > oo limit. A spacelike geodesic with a = 1 will approach a boundary null 
geodesic in the limit. Its end point can be anywhere in the straight line between the origin 
and (it, it) depending the value of E — J which should be fixed in the limit E,J—> oo. 
A spacelike geodesic with a G (l,ao) will approach a bulk null geodesic in the limit. A 
spacelike geodesic with a > ocq will not come back to the same boundary. Fig. 21 plots the 
end points of spacelike geodesies with E/J fixed. 



Appendix D: No coming back for Poincare patch 

In the following we argue that for field theories formulated on R^" 1 ' 1 and states respecting 
the full Poincare symmetry there are no new light-cone singularities. As discussed in Section 
2, the presence of light-cone singularities in the boundary field theory is governed by the 
properties of bulk null geodesies. For field theories formulated on R d_1,1 we will show that 
there are no null geodesies through the bulk spacetime connecting boundary points. This 
implies that the only null geodesies that connect points on the boundary are those that lie 
entirely within the boundary and hence the nature of singularities in the correlation function 
is entirely determined by the boundary causal structure. 

To establish the absence of null geodesies connecting boundary points through the bulk, 
let us consider a bulk spacetime with negative cosmological constant foliated by R d_1,1 slices. 
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The metric can be written in the warped-product form by picking a metric g m n{x l ) on the 
boundary and a radial coordinate r: 

ds 2 = e 2A(r) g mn dx m dx n + dr 2 (D. 1) 

For this geometry one can write down the geodesic equations as 

•2 , 2A(r) ~m-n _ 

I I O- L/ ,j'f'i yi Jii JU (V . 

(D.2) 

r-A'(r)e 2A ^g mn x m x n = 

with k being ±, for spacelike, timelike and null geodesies, respectively. By eliminating the 
boundary directions we can write an effective classical particle in a potential equation for 
the motion in the radial direction: 

r + A'(r) (r 2 - «) = (D.3) 

The issue of whether geodesies emanating from the boundary into the bulk turn around, 
can now be analyzed simply in the effective problem for the radial motion (D.3). For the 
case of null geodesies which we are especially interested in the first integral from (D.2) gives: 

(D.4) 



D A(r) 



implying that A(r) —* oo for the geodesic to turn around (since r = by definition at the 
turn around point). 

However, A(r) cannot diverge to oo - the energy condition on the stress tensor require 
that A'(r) be a monotonically decreasing function [32]. This effectively implies that A(r) 
decreases from oo at the AdS boundary where the field theory is formulated. Physically 
divergence of A(r) would look like we have second timelike boundary in the geometry which 
again is disallowed by the null geodesic convergence condition in spacetimes where the matter 
sources respect the null energy condition. 

Note that this result is contingent on the Poincare symmetries of the boundary being 
preserved. 30 We can have geodesies that return to the boundary even in the non-compact case 
when we consider states that break the boundary Poincare invariance. A simple example is 
to consider the global AdS-Schwarzschild geometry that is sliced in the Poincare coordinates 
(c/., [1] for a boundary description of this state). Here it is possible to exploit the black 
hole, which breaks the Poincare symmetry to slingshot oneself back to the boundary in a 
timescale shorter than the AdS time. 



30 We thank Simon Ross for pointing this out to us and suggesting the example described. 



42 



Appendix E: Collapsing shell in AdS 



In Section 5 we derived a scaling formula for t Q as a function of 5t = th — U using a 
ray tracing calculation. The exponential relation is given in (5.6) and intuitively we can 
imagine this arising from the red-shift in the vicinity of the horizon. We now proceed 
to derive this result more explicitly, by finding actual radial null geodesies in the collapse 
background. In this Appendix we constrain ourselves to the sharp shell cleanly separating 
AdS and Schwarzschild-AdS; in Appendix F we consider the more general smeared shell 
collapse modeled by a Vaidya spacetime and also extend our analysis to non-radial geodesies. 

Consider the set-up as sketched in Fig. 8. There are two interesting limits to consider 
as discussed in Section 5.3. In particular, we wish to obtain the leading behaviour both as 
ti — > t^ and as tj — > t s , as indicated in Fig. 12a and b, respectively. In fact, in order to 
derive both (5.6) and the scaling exponent 7 introduced in Section 5.3, we will find the full 
exact expression for t a in terms of ti and the parameters t s and r + describing the shell. It is 
interesting to examine the geodesic behavior both in d = 3 and in d = 5. The distinction of 
course is engendered by the fact that the BTZ black hole in d = 3 being an orbifold of AdS3 
is a simpler geometry. On the other hand, the physics of the type III geodesies of Fig. 8b is 
different is quite different since spacelike geodesies do not bounce off the BTZ singularity, 
so we can meaningfully consider only the ti — > scaling. 

The strategy for calculating radial null geodesies will be to write the geodesies in the 
different coordinate patches (such that in each patch we can use the metric of the form in 
(5.1)), and then patch them together using the fact that both the geodesies and the shell 
are null. The three regions of interest are depicted in Fig. 12c, and correspond to 

1. pure AdS (before/inside the shell) 

2. Schwarzschild-AdS, outside the horizon (outside the shell) 

3. Schwarzschild-AdS, inside the horizon (but still outside the shell) 

E.l d = 3 

As a warm-up, let us first focus on d = 3. For simplicity, we use the metric (5.1) with 
f(r) = r 2 + 1 in region 1, and f(r) = r 2 — r 2 + in regions 2 and 3. We will write the 
expressions for geodesies in 3 distinct regions (even though in 3 dimensions, region 3 is 
not so relevant since radial spacelike geodesies in BTZ do not bounce off the singularity). 
Calculating the geodesies as before, we have for pure AdS: 



tAds(r) =t ± (tan 1 r - tan 1 r ) 



(E.l) 



and for BTZ: 




inside 
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where we denote the initial conditions by t(r ) = t - We can now patch these together to 
see when a geodesic starting at t = ti, r = oo reemmerges back out to r = oo. Denoting this 
time t Q , we can express t Q in terms of ti, t s , and r+, obtaining: 

1 /tan(^ i )+r + \ 
'^'^^'" Ut^-rJ (E3) 

Matching where the shell intersects the horizon, we can relate t s to th as follows: 

t s =t h + 2 tan" 1 r+ (E.4) 

Finally, substituting (E.4) into (E.3) and expanding for small 5t, we obtain the scaling 
behaviour 

t -In \8t\ 5t e~ r+t ° = -e~ Kt ° as t -> oo (E.5) 

r + 

where k = r + is the surface gravity of the BTZ black hole. 

Unfortunately, in this set-up we can only answer the first part of the questions raised 
in Section 5.3: as 5t — > 0~, t Q diverges logarithmically with the coefficient given by the 
black hole temperature, as one would naively expect. To see the other scalings, t; t — > t\ and 
U — > t s , discussed in Section 5.3 , we need to go to higher dimensions. 31 



•+r + J \r 



outside 



(E.2) 



E.2 d = 5 

Now let us consider d = 5. Whereas the time-radius relation for radial geodesies is indepen- 
dent of the dimension in pure AdS, we have a more complicated relation t(r) for the black 
hole geometry than in the BTZ case. Writing 

t(r) = t(r ) ± jf (E.6) 

and rewriting the metric function outside the shell as 

/(r) = (r 2 + ^ + )(r 2 -rj) (E ?) 

where p 2 + = 1 + r 2 + , we can express the indefinite integral pertaining to outside and inside 
horizon regions respectively as 



f dr If i r i r + \ 

/ = p + tan r + tann — for r > r + 

J f{r) r + n V P+ r / 



;e.s) 



31 If we nevertheless ignore the fact that spacclike geodesies do not bounce off the singularity, and consider 
the two null geodesies meeting at the singularity, then we would expect to find 1) the same scaling behaviour 
(same coefficient) as in (E.5) for ti > th, and 2) (t s — ti) ~ (t a — t s ) (i.e. 7 = 1). 
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/dr If , r r \ 

= p + tan r + tanh — tor r < r + 

/(r) r + K\ p+ r + J 



(E.9) 



X_i_2r 2 

where k = r+ + is the surface gravity of the black hole. For a geodesic crossing the horizon, 
the infinite contribution from the coordinate singularity cancels out, so if T\ > r + and 
r 2 < r + , say, we would obtain 



tsAdsfa) = t S Ads(ri) — ( p+ tan 1 — 



T2 

+ r+ tanh^ 1 — - r+ tanh -1 — ] (E.10) 



n 



Applying this to the geodesies of type II and III in Fig. 12a, we obtain the full exact 
expression for t in terms of t s , the parameters specifying the black hole, and the radius of 
where the shell crosses our geodesic r c : 

2 / 7T T V \ 

t a = t s -\ I p + p + tan -1 — + r + tanh - — ) for geod. II (E.H) 

r + K V 2 p + r + J 

2 / r r \ 

t a = t B -\ -p+ tan -1 — + r+ tanh -1 — for geod. Ill (E.12) 

r + K \ p + r + J 

where (in both cases) the crossing radius is given by 

r c = tan r^i\ (E.13) 

and the time of the shell is related to the time of horizon formation by 

t a = t h + 2 tan" 1 r+ (E.14) 

We can easily expand this out for small 5t, to check that in both cases 

\5t\xi e~ Kt ° t ~-l nr — . (E.15) 

K \0t\ 

This is again consistent with our expectations (5.6), and indeed of the same form as the 
result obtained for the 3-dimensional case. 

Finally, we can also use expression (E.12) to consider the case of Fig. 12b. In particular, 
if we let t s —ti = e, then expanding t — t s to third order in e, we find that 

to-ts^^-^e 3 . (E.16) 

iz / + p + 

This in particular implies that the scaling parameter 7 introduced in Section 5.3 is 3 for 
black hole formation in d = 5. Intriguingly, we can generalize this to d dimentions: For 
<i-dimensional collapse spacetime, 7 = d — 2. 
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Appendix F: Smeared shell in AdS and non-radial geodesies 

In Appendix E, we have considered radial null geodesies in the collapsing null shell 
geometry. In particular, we calculated t Q as a function of U and the shell parameters. As 
discussed in Section 5, being able to read off the value of ti for which t Q diverges from the 
gauge theory correlators, we can automatically determine the horizon formation time th 
(=>- t H ). Here we wish to extend this analysis in two directions. First, we want to consider 
a more general (and physically more realistic) background corresponding to a smeared shell. 
Second, we want to analyze non-radial geodesies, as motivated in Section 5. To that end, we 
will discuss general null geodesies in Vaidya-AdS spacetime. 



F.l Vaidya for AdS 

Consider the stress tensor for a spherical null gas 

T w = h(r,v) (F.l) 

The metric can be written as 

ds 2 = -f(r,v)dv 2 + 2dvdr + r 2 dn 2 3 (F.2) 

with 

/M) = r 2 + 1-^ (F.3) 
The covariant conservation V M T Ml/ = of T implies that 32 

h(r, v ) = — q(v) (FA) 
The Einstein equations now reduce to (with A = —6) 



2 

m'(v) = —q(v) (F.5) 
For the sharp shell we take q(y) has the form of a delta function, i.e., 

m(y) = v < Vq (F-6) 
= n v>v (F.7) 

with q{y) = ^fi5(v — v ). Note that in this case the w-dependence disappears before and 
after the shell; the spacetime (F.2) is static in both regions (but not globally static because 
of the shell), written in ingoing coordinates, 

dv 

v = t + z, dz = Y 



:;2 



Note that T v rv = = 0. 
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F.2 Null geodesies 



Let us now consider null geodesies in the background (F.2). Writing the second-order 
geodesic equation (in the equatorial plane) directly yields 

v + \ d r f(r, v) v 2 — r ip 2 = 
r - \ [d v f(r, v) + f(r, v) d r f(r, v)} v 2 - r f(r, v) (p 2 = 

</3 + fr0 = O (F.9) 

For any given f(r, v), we can solve these numerically to find any geodesic through the bulk. 
The function f(r, v) with which we choose to model the smeared shell in Section 5 is 

a fl + tanh^\ 
f(r, v) = r 2 + 1 - £ ( ^\ (F.10) 

which smoothly interpolates between AdS at v — > — oo and Schwarzschild-AdS as v — >• +oo. 
In particular, the shell is inserted at 33 v = = t s + ir/2, and has 'thickness' v s . This provides 
a convenient regulator of the thin shell collapse, since as v s — >• 0, we recover the collapse 
spacetime (5.1) written in ingoing coordinates. 

However, to find the geodesies in the collapse spacetime, and specifically to determine 
th{o)i we can use a simpler and more elegant method, which we now indicate. The geodesies 
can be obtained from the following action 

S = J dr {- l -fv 2 + vf + irV) (F.ll) 

The canonical momenta are 

J = r 2 <P, E = -°— = -r + fv, °— = v (F.12) 
ov or 



Note that J is conserved, while E is in general not. Note that J can be scaled to be 1 by a 

E 
J 



rescaling of r. As before, we will denote a = ^ and set J = 1. Another first integral of the 



system is given by 



r 2 + 4 = « 2 (F-13) 



r 



The equations of motion can now be written as 

19/ ,2 



a = (F.14) 



33 Since both AdS and Schwarzschild-AdS geometries are static, we can WLOG set the time t s of the shell's 
creation on the boundary; for convenience we choose t s = —tt/2. 
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From the above equations we can also derive 



where we have used that 



ldf 2 1 df f . i 

f = — -v 2 + — (F.16) 

2 dv 2r 2 dr r 3 



a + r 

(F.17) 



/ 



Now we apply the above equations to the Vaidya spacetime with a sharp shell, in which 

case 

^ = -^<%— o) (F-18) 

The discontinuity of / across v = v is thus given by 5f = —^3. One can readily conclude 
from the above equations that both a and f jump across v = v , in fact by the same amount, 
while v is continuous. More explicitly we find that 



5r — 5a — — - -7 „ u 



< (F.19) 



Let us now return to the problem of determining th(a) for the collapse geometry. Re- 
alizing that a jumps across the shell, let us work backwards by considering what feature 
of the geodesies makes t large. The relevant spacetime to consider is Schwarzschild-AdS, 
wherein the geodesies follow from the potential (4.5) drawn in Fig. 6. Since the top of 
the potential is always outside the horizon (as follows from the discussion in Section 4, 

= 2/i = 2r 2 + (7 + + 1) > r+), the only way that a non-radial null geodesic could emerge at 
infinite time would be to get trapped in the unstable circular orbit at the top of the effective 
potential. This requires 

a 2 = V eff (r m ) = l + ±- (F.20) 

Knowing the a in the Schwarzschild-AdS part of the spacetime, we can now use the expression 
(F.19) for the jump in a across the shell, to find what initial conditions (i.e. what in AdS) 
we need to start with to achieve this a . Using (F.17), (F.13), and the expression for f(r) in 
AdS (5.2), all evaluated at some crossing radius r x where the geodesic intersects the shell, 
we can reexpress (F.19) as follows: 



a ° = a *-2^ x — — (R21) 

Finally, we can express the crossing radius r x in terms of the initial time by following the 
geodesic with parameter a.{ from r x back to its starting point on the boundary at tj. This is 
given by 

t s -U = tan" 1 r x + tan" 1 l - — (F.22) 

Oti 
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Note that ti should be viewed as a function of ati, and for t Q — > oo, this is nothing but th{a) 
of Section 5, where we have been using a = «j. 

Hence, to find th{ct), we can solve the system of equations (F.20), (F.21), and (F.22). 
We first solve (F.20) and (F.21) for the crossing radius r x (ai), and then substitute this into 
(F.22) to find th(ce) = ti, as a function of a, = a and the shell parameters, \x and t s . This 
determines the black curve plotted in Fig. 9. Note that the radial geodesic limit (a — > oo) 
is continuous, despite the fact that only the radial geodesic can truly sample the horizon 
formation event, since only the null outgoing radial geodesies can escape from the close 
vicinity of the horizon. 
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